Spin-dependent resonant tunneling in symmetrical double-barrier structure 
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A theory of resonant spin-dependent tunneling has been developed for symmetrical double-barrier 
structures grown of non-centrosymmetrical semiconductors. The dependence of the tunneling trans- 
parency on the spin orientation and the wave vector of electrons leads to (i) spin polarization of 
the transmitted carriers in an in-plane electric field, (ii) generation of an in-plane electric current 
under tunneling of spin-polarized carriers. These effects originated from spin-orbit coupling-induced 
splitting of the resonant level have been considered for double-barrier tunneling structures. 

PACS numbers: 72.25.Dc, 72.25.Mk, 72.25.Hg, 73.40.Gk 
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I. INTRODUCTION 



Physics of spin-dependent tunneling phenomena in 
semiconductor structures has attracted lately a great deal 
of attention. Significant progress has been made in ex- 
perimental and theoretical investigation of spin-polarized 
transport in magnetic tunneling junctions (for review 
see [IB]). On the other hand, it was pointed out re- 
cently that electron tunneling could be spin-dependent 
even in the case of nonmagnetic barriers. It was demon- 
strated that the transparency of a semiconductor barrier 
depends on the spin orientation of carriers if the sys- 
tem lacks a center of inversion. Two microscopic mecha- 
nisms, Rashba spin-orbit coupling induced by the barrier 
asymmetry^*^ and the fc 3 Dresselhaus spin-orbit split- 
ting in non-centrosymmetrical materials?^ were shown 
to be responsible for the effect of spin-dependent tunnel- 
ing. Spin-orbit interaction couples spin states and space 
motion of conduction electrons that opens a possibility to 
orient, manipulate and detect spins by electrical means. 
Effect of spin-dependent tunneling in nonmagnetic semi- 
conductor heterostructures was supposed to be applied 
for spin inject ior&2*i2*ii and detection of spin-polarized 
carriers^li Devices based on spin-dependent tunneling 
were suggested to be utilized as components of the spin 
field-effect transistor^ 



In this paper we present a theory of spin-dependent 
tunneling through a double-barrier structure, the struc- 
ture of resonant a tunnel diode (RTD), grown of non- 
centrosymmetrical semiconductors. Section^is devoted 
to calculation of the spin-dependent transmission coeffi- 
cient. The dependence of the structure transparency on 
the spin orientation and the wave vector of electrons can 
be employed for spin injection and detection: (i) an elec- 
tric current flow in the plane of interfaces leads to the 
spin polarization of the transmitted carriers, (ii) trans- 
mission of the spin-polarized carriers is accompanied by 
generation of an in-plane electric current. These effects 
are considered in Sections IIIII and llVl respectively. The 
results of the numerical calculations are compared with 
that obtained in a simple analytical theory. 



II. SPIN-DEPENDENT RESONANT 
TUNNELING 

We consider the transmission of electrons with the ini- 
tial wave vector k = (fey , k z ) through a symmetrical dou- 
ble barrier structure grown along z || [001] direction (see 
FigQJ. Here fen is the wave vector in the plane of the in- 
terfaces, and k z is wave vector component normal to the 
barrier and pointing in the direction of tunnelling. The 
electron motion in each layer of the structure is described 
by the Hamiltonian 



H 



h 2 k\ 

2m* dz 2 2m* 



d 2 



+ V(z) + H D , 



(1) 



where m* is the electron effective mass, V(z) is the het- 
erostructure potential equal to Vb for the barriers and 
—V w for the well, and Tiu is the spin-dependent k 3 Dres- 
selhaus term that describes spin-orbit splitting of the con- 
duction band in zinc-blende-lattice semiconductors. We 
assume the RTD structure to be designed so that the res- 
onant transmission occurs for the incident electrons with 
the kinetic energy e much smaller than Vb and V w . Then 
the Dresselhaus term is simplified tc 7 
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(2) 



where 7 is a material constant, a x and a v are the Pauli 
matrices, and the coordinate axes x,y,z are assumed to 
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FIG. 1: Electron transmission through (OOl)-grown RTD 
structure. Vb is the height of the barriers, V w is the depth 
of the well, a and b are the well width and barrier thickness, 
respectively. 
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be parallel to the cubic crystallographic axes [100], [010], 
[001], respectively. 
The Dresselhaus term (0 is diagonalized by the spinors 



X± 



~( 1 

V2\ Te- 



rn 



which describe the electron states "+" and "— " of the 
opposite spin directions. Here ip is the polar angle of the 
wave vector k in the xy plane, fc|| = (fc|| cos<^ , k\\ sirup). 
The electron spins s± corresponding to the eigen-states 
"±" are given by 

s±{k\\) = ^X±^X± = \ {Tcoap, ±siny>, 0) . (4) 

One can note that essentially Hp induces a spin- 
dependent correction to the effective mass for electron 
motion along z. In the basis of the spin eigen-states "±" 
the effective Hamiltonian has the simple form 



H± = 



d 2 



2 1.2 



2m± dz 2 2m* 



V(z), 



(5) 



where the effective mass along z, modified by spin-orbit 
coupling, depends on the in-plane electron wave vector 
and is given by 



m± = m 



1±2 



jm*k\\ 



(6) 



Solution of the Schrodinger equation with the Hamilto- 
nian JSJ and boundary conditions for the wave functions 
■0±i which require the continuity of 



and 



1 dip± 
m± d z 



at the interfaces, allows one to derive coefficients of trans- 
mission, t±, and reflection, r±, for the electrons of spin 
eigen-states "+" and "— ". 

FigEl presents the dependencies of the double-barrier 
structure transparency \t±\ 2 on the incident electron en- 
ergy along growth direction, e z = h 2 k 2 /2m*, calculated 
numerically for the fixed in-plane wave vector fen . The 
spin splitting of the resonant peak is clearly seen. In 
calculations both the electron effective mass m* and the 
Dresselhaus constant 7 are assumed to be the same for 
the barrier and well layers. 

In the limit of thick barriers the structure transparency 
demonstrates sharp peaks and hence can be approxi- 
mated by Dirac (^-functions 

\t±(e z ,h)\ 2 « nT ± (k {l )d[e z -E ± (k {[ )], (7) 

where the prefactors T±(k^) describe the transmission 
efficiencies, and E±(k\\) stand for the energies of the res- 
onances. The positions of the resonances E±(k\\) corre- 
spond to the energies of size-quantization of an electron 
in the quantum well of infinitely thick barriers with the 
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FIG. 2: The transparency of the double-barrier structure, 
jt±| 2 , as a function of e z at fixed fcii = 10 6 cm -1 . The in- 
sert shows the dependence of the spin splitting of the res- 
onant peak on ku calculated numerically (solid curves) and 
following Eq.® (dashed lines). The used parameters, 7 = 
76 eV-A 3 , m* = 0.053m , V b = 230 meV, V w = 200 meV, 
a = 30A, and b = 50A, correspond to AbGai-^Sb, x = 
0.15/0.3/0/0.3/0.15, RTD structured 



Dresselhaus spin-orbit interaction included. The trans- 
mission efficiencies r±(fc||) are determined, to the con- 
trary, by the electron lifetimes on the resonant levels in 
the double-barrier structure. Considering the spin-orbit 
interaction to be small perturbation, the positions of the 
resonant levels and their widths can be expanded as 

E ± (k l{ ) = E ± ak l{ , r±(fc||) = (l±/3ft]|)r , (8) 

where Eq and Tq are the level position and the width 
when the spin-orbit interaction is neglected, 



T = 8y/E Q (V b -E Q 



(V b - E )(V W + E ) exp(-2 K &) 
V b (V b + V W ) 1 + na/2 



(9) 

k = y/2m*(Vb — E )/h is the reciprocal length of the 
wave function decay under the barrier, a is the quantum 
well width, b is the thickness of the barriers, and a and 
(3 are the coefficients which describe the spin splitting of 
the level and the width modification of the spin sublevels, 
respectively. For the case under study, Eq <g; Vb,V w , they 
are given by 



27m* 



1 + 2/na 



(10) 
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Note, that the width To and the coefficient (3 are highly 
sensitive to the position of the resonant level Eq . For the 
RTD structure parameters presented in the caption to 
FigEl the coefficients can be estimated as follows: a = 
9.7-10- 7 meV-cm, j3 = (4.2-10~ 8 +4.8-10^ 7 -meV/i;o) cm. 
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III. SPIN ORIENTATION OF CARRIERS 

Spin splitting of the resonant peak at non-zero feii can 
be employed for injection of spin-polarized carriers. We 
assume two parts of bulk semiconductor separated by the 
RTD structure, and electrons tunneling through the bar- 
rier from the left to the right side of the structure. In 
equilibrium the momentum distribution of the incident 
electrons is isotropic in the interface plane and there- 
fore the average spin of the transmitted carriers vanishes. 
This isotropy can be broken by application, for example, 
of an in-plane electric field, F . Then the carriers tunnel 
through the structure with non-zero average wave vector 
in the plane of interfaces that leads to the spin polariza- 
tion of the transmitted electrons*^ 

The average spin of the transmitted electrons is given 

by 



s = S/N , 



(11) 



where S and N are the spin and carrier fluxes through 
the barrier given in the linear in F regime by 



s= /i(fe) fen)| 2 s + (fe||) 



(12) 



fell , kz>0 



+ |i-(e 2 ,fc||)| s-(fey) 



|t+(ea,fe||)| + |t_(e a ,fc||)| 



fen , k z >0 



v = hk/m* is the electron velocity, /o(e) is the equilib- 
rium distribution function, e is the electron energy, fi(k) 
is the electric field-induced correction to the distribution 
function, 



-er p — („„ • F) 



e is the electron charge, and r p is the momentum relax- 
ation time. 

Substituting the spin vectors of the Dresselhaus eigen- 
states s± in the form Eq. J3J into Eq. I)12[l. one derives 
the angular dependence of the average spin of the trans- 
mitted carriers, 



Vd,: 

2v d 



P, 



(13) 



where Vd — (eT p /m*)F is the in-plane drift velocity of 
the incident electrons, and P s is the spin polarization of 
the transmitted particles, 



P, 



Vdm 
2N 



E 

fen , k z >0 



dfo 
de 



|*+( £ *)fc||)| ~ |*-( £ 2! fc ||)| 



V Z V\\ 

(14) 

The direction of the electron spin (|13|l is determined by 
the symmetry of the Dresselhaus term. In particular, the 
spin s is parallel (or antiparallel) to the electron drift 



velocity Vd, if I'd is directed along the crystal cubic axis 
[100] or [010]; and s is perpendicular to v d , if the latter 
is directed along the the axis [110] or [110]. 

In the limit of thick barriers the structure transparency 
can be approximated by Dirac 8- functions. Then substi- 
tuting the transparency |t±| 2 in the form Eq. J7J and 
assuming spin-orbit interaction to be small, one derives 
the following expression for the spin polarization of the 
transmitted carriers 



u VdTn i it as 



(15) 



where C = fo(s)de / fo(E ) is an energy equal to 
Ep — Eq for 3D Fermi and fc^T for 3D Boltzmann gas, 
respectively. 

The dependence of the spin polarization P s on the res- 
onant level position Eq is presented in FigOl fo r the in- 
cident electrons forming degenerate gas with the Fermi 
energy Ep = 10 meV and the gas of the same carrier 
concentration at temperature T = 300 K. Solid curves 
correspond to the exact solution of Ea. p4|) . with the 
transmission coefficients t±(s z ,h\) being calculated nu- 
merically in the framework of the transfer matrices. The 
variation of Eq from to 10 meV was achieved modify- 
ing the quantum well width from a « 3lA to a w 29.4A. 
Dashed curves are plotted following Ea.pSjl for Fermi 
and Boltzmann statistics. One can see that in the wide 
range of the resonant level positions the simple analytical 
theory demonstrates a good agreement with the numer- 
ical calculations. For the reasonable set of parameters 
given in the figure caption it is possible to achieve spin 
polarizations of several percents. The sign of the po- 
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FIG. 3: The spin polarization P 3 as a function of the reso- 
nant level position Eq for degenerate electron gas with Ef = 
10 meV and non-degenerate gas of the same carrier concentra- 
tion at T — 300 K. Solid curves correspond to the numerical 
calculation, dashed curves are plotted following Ea. l|15^ . The 
parameters of the double-barrier structure are presented in 
the caption to Fig|5| and vj = 2.5 ■ 10 6 cm/s that corresponds 
to Vf/W. 
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larization P s is governed by interplay between the spin 
splitting, afcy , and the difference of the widths of the spin 
sublevels, /3fc|| , since the carrier occupation of the lower 
sublevel "— " is larger than that of the higher sublevel 
"+" while the tunneling transparency of the sublevel "— " 
is smaller than that of "+" (see Eq.lJJJ)). It clarifies why 
the terms proportional to a and /? contribute to Eq. l|15|) 
with opposite signs. At low temperature both terms are 
comparable, and P s changes the sign with increasing of 
E . At E <C Ep the effect is mainly determined by the 
difference of the spin sublevel transparencies. The role 
of the energy spacing between the sublevels is negligi- 
ble, since the carrier populations at the energies E+ and 
E- almost coincide. With increasing of Eg the "step- 
like" Fermi distribution leads to different occupations of 
the spin sublevels, and the polarization P s changes the 
sign. At higher temperatures the carrier distribution be- 
comes to be smooth, and the effect is related mainly to 
the difference of the tunneling transparencies of the spin 
sublevels. 



where fo is the equilibrium distribution function of non- 
polarized carriers, n s is the unit vector directed along 
the spin orientation, p s is the degree of the polarization, 
and (1/e) is the average value of the reciprocal kinetic 
energy of the carriers equal to 3/Ep for 3D degenerate 
electron gas with the Fermi energy Ep, and 2/ksT and 
3D non-degenerate gas at the temperature T. The spin 
matrix of the electron transmission through the structure 
is given by 



T = t + (e z , % ) x+X+ + t- (e z ,h) X -X- 



(18) 



Substituting the density matrix l|17|l and the transmis- 
sion matrix (|18(l into Ea. (|16|) and taking into account 
the definition of the vectors s± J2J, one obtains 



4er p p s 



17 |M £ *> fc l|)| 2n «' s +( fc l|) 



(1/e) ^ de 



,h)\ n s 



S_(fc|| 



(19) 



IV. TUNNELING SPIN GALVANIC EFFECT 

Generation of an electric current by spin-polarized car- 
riers represents the effect inverse to spin injection. Now 
we assume the electron gas on the left side of the double- 
barrier structure to be spin-polarized. Electrons with 
various wave vectors tunnel through the RTD. However, 
due to the splitting of the resonant level the flux of the 
spin-polarized carriers with the certain in-plane wave vec- 
tor fen is larger than the flux of the particles with the 
opposite in-plane wave vector, — fey. This asymmetry re- 
sults in the in-plane flow of the transmitted electrons near 
the barrier, i.e. in the interface electric current. The di- 
rection of this interface current is determined by the spin 
orientation of the electrons and symmetry properties of 
the barrier, in particular the current reverses its direction 
if the spin orientation changes the signpi2iH 

The theory of such "tunneling spin-galvanic effect" is 
developed by using the spin density matrix technique. 
The interface current of spin-polarized electrons trans- 
mitted through the tunneling structure is given by^ 



3\\= e ^2 T p v \\ v z^r[TpiT^] 



(16) 



fe||,fcz>0 



where r p is the momentum relaxation time, pi is the elec- 
tron density matrix on the left side of the structure, and 
T is the spin matrix of the tunneling transmission that 
links the incident spinor wave function ipi to the trans- 
mitted spinor wave function ip r , ip r = Tipi. We assume 
the carriers on the left side of the structure to form 3D 
spin-oriented electron gas, and electron distributions in 
both spin subband to be thermalized. For the case of 
small degree of spin polarization, the density matrix has 
the form 

dfo 2p s 



Pi « fol 



de (1/e) 



(17) 



Taking into account the explicit form of the spin vec- 
tors of the Dresselhaus eigen-states (@J), the components 
of the tunneling spin-galvanic current have the form 



J\\,x = ~J\\ 



31 



JW 



(20) 



where 



n 



eTpp s \ - 
"(1/e) ^ 



dfo 
de 



VzVu 



The direction of the tunneling spin-galvanic current is 
determined by the spin orientation of the electrons with 
respect to the crystallographic axes. 

In the limit of thick barriers the structure transparency 
can be obtained substituting (5-functions (JJJ for the struc- 
ture transparency. Then assuming spin-orbit interaction 
to be small, the interface current is derived to be 



n 



eT p p s f(E )m* 
(1/e) nh 4 



(22) 



Fig El presents the dependence of the interface current 
ju on the energy position of the resonant level E for the 
incident electrons forming spin-polarized degenerate gas 
with the Fermi energy Ep = 10 mcV and the gas of the 
same carrier concentration at temperature T — 300 K. 
Solid curves correspond to the exact solution of Eq. 1 12 11 1 , 
with the transmission coefficients t±(e z ,k\\) being cal- 
culated numerically. Dashed curves are plotted follow- 
ing Ea. l22|) for Fermi and Boltzmann statistics. For 
the AlGaSb-based RTD structure considered as an ex- 
ample the tunneling spin-galvanic current is of order of 
mA/cm. Estimation shows that it is enhanced by an 
order of magnitude with respect to the interface cur- 
rent generated under tunneling through a single AlGaSb- 
based barrier provided the equal electron tunneling flux, 
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FIG. 4: The interface current jn as a function of the reso- 
nant level position Eo for degenerate electron gas with Ef ~ 
10 meV and non-degenerate gas of the same carrier concentra- 
tion at T — 300 K. Solid curves correspond to the numerical 
calculation, dashed curves are plotted following Eq.!^. The 
parameters of the double-barrier structure are presented in 
the caption to Fig[2] t p = 1 ps, and p a = 0.1. 



N ~ 10 22 l/(cm 2 s)ii^ Similarly to the spin polariza- 
tion (Fig[3J), the sign of the interface current (1221) is gov- 
erned by interplay between the contributions responsible 
for the spin splitting and the difference of the spin sub- 
level transparencies. 



In conclusion, the theory of spin-dependent electron 
tunneling has been developed for symmetrical double- 
barrier structures based on zinc-blende-lattice semicon- 
ductor compounds. The Dresselhaus spin-orbit interac- 
tion couples spin states and space motion of conduction 
electrons that leads to spin splitting of the resonant level 
depending on the in-plane electron wave vector. The ef- 
fect of the spin-dependent tunneling could be employed 
for creating spin injectors and detectors based on non- 
magnetic tunneling structures. 
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